Universal Conductance and Conductivity at Critical Points in Integer Quantum Hall Systems 
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The sample averaged longitudinal two-terminal conductance and the respective Kubo-conductivity are cal- 
culated at quantum critical points in the integer quantum Hall regime. In the limit of large system size, 
both transport quantities are found to be the same within numerical uncertainty in the lowest Landau band, 
0.60 ± 0.02 /h and 0.58 ± 0.03 e^/Zi, respectively. In the 2nd lowest Landau band, a critical conductance 
0.61 ± 0.03 /h is obtained which indeed supports the notion of universality. However, these numbers are 
significantly at variance with the hitherto commonly believed value 1/2 e^//i. We argue that this difference 
is due to the multifractal structure of critical wavefunctions, a property that should generically show up in the 
conductance at quantum critical points. 



I 

O 

o 



> 
o 



in 
o 



I 

O 

o 



PACS numbers: 72.15.Rn, 71.30.+h 

The transitions between quantized Hall plateaus are gener- 
ally viewed as a manifestation of continuous T = quantum 
phase transitions. They emerge in the presence of a strong 
magnetic field in disordered two-dimensional electron or hole 
gases which can be observed in Silicon MOSFETs and GaAs- 
heterostructures. According to renormalization group consid- 
erations, some properties of the quantum phase transitions 
should be universal, i.e., independent of microscopic details 
of the system. For instance, this is supposed for the criti- 
cal exponent that controls the divergence of the localization 
length £,{E) = £,o\E — E^^^ near the critical energies E^^a 
behavior which has been corroborated by experimental []J,|2|] 
and numerical jsl 0, 01 studies. In addition, it has been sug- 
gested that the impurity averaged elements of the conductivity 
tensor act as coupling constants in an appropriate field theory 
0] where all the transitions belong to the same universality 
class 0] with universal values of the critical dissipative and 
Hall conductivities. It has been proposed |8, 9] that their val- 
ues (j^y = (n — 1/2) e^/h and cr^^ = l/2e^//i are simply 
related at the respective 7i-th transition. 

Away from the critical points, the Hall conductivity is ex- 
actly quantized due to topological reasons. Thus, in the 
presence of electron-hole symmetry, the above formula for 
axy is evident at the transition points fioll . Support for 
the existence of a universal axx comes both from experi- 
ments 

El El and numerical calculations. Amongst the 
latter, a critical conductivity a^^ ~ 0.54 e^/h was found 
for the lowest Landau band using white noise disorder po- 
tentials E3I1 . A similar diagonalization study Eoll yielded 
'^xx — (0-55 ± 0.05) e'^/h for various correlated random po- 
tentials, a result considered by the authors to be consistent 
with the prediction 0] a^^ — 1/2 e'^/h. However, the critical 
diagonal conductivity deduced from a Thouless number study 
E4] exhibited a strong dependence on the range of the scat- 
terers. Moreover, calculations of the two-terminal conduc- 
tance of square systems produced a significantly larger value 
gc — (0.58± 0.03) e^/ h |, I.S | and showed considerable fluctu- 
ations exhibiting a broad distribution between zero and e'^ /h 
E3l. Here gc = lim^^oo J2i 9ii^)/-^ is the disorder aver- 



aged critical conductance of N realizations. This value was 
derived from a finite-size scaling study of the transfer-matrix 
results calculated within a Chalker-Coddington (CC) network 
model Q]. In a recent work, this outcome has been confirmed 
with gc = (0.57 ± 0.02) e^/h El- 

On the other hand, similar transfer-matrix calculations 
within a tight-binding (TB) lattice model with spatially uncor- 
rected disorder potentials yielded gc = (0.506 ± 0.01) e'^/h 
E^ . hence, the asserted universality cannot be taken for 
granted. Even if it were possible to reconcile the contradict- 
ing results, the hitherto unknown relation between the two- 
terminal conductance and the conductivity, the latter calcu- 
lated for instance from the Kubo formula for systems without 
contacts and leads, would still inhibit us drawing final conclu- 
sions about the validity of the proposed universal value of the 
critical dissipative conductivity tj^^ = 1/2 e'^/h H. 

Our work is aimed at resolving this issue by making avail- 
able the missing links and therewith provide evidence for 
an equality of conductivity and conductance at the critical 
point in the lowest Landau band. We find, within the un- 
certainty of our calculations, similar values 0.58 e^/h and 
0.6e^/h, respectively. These values are, however, signifi- 
cantly larger than the previously proposed and well accepted 
a^x " l/2e^//i. The origin of this discrepancy is suggested 
to arise from the multifractal correlations of the critical eigen- 
state amplitudes, which were presumably not fully taken into 
account previously. We also present the critical conductance 
distribution for the second Landau band which agrees with the 
one from the lowest. The scaling of the ensemble averaged gc 
yields (0.61 ± 0.03) e'^/h, a result that supplies further evi- 
dence for the existence of a universal conductance. 

We first show that for non-interacting electrons moving in a 
quenched disorder potential and perpendicular magnetic field 
described by a tight-binding lattice model the following rela- 
tion holds. 



gc{L) = 5c(oo) - go {Lo/L)^, 



(1) 



if the size of the square systems is L > Lq. We find in the 
lowest Landau band a critical conductance gdoo) = (0.60 ± 
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0.02) e^/h with y = 0.4 ± 0.02. Such a relation has been 
suggested previously fl9ll where the least irrelevant scaling 
exponent y was proposed to equal the exponent of the critical 
eigenfunction correlations 11311 . Then, we demonstrate that 
the Kubo-conductivity for very long strips of width M > Mq 
obeys a similar power law 
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,(oo)-ao(Mo/M)^ 



(2) 



with (^^^.(cx)) = (0.58±0.02) e^//i. The irrelevant scaling ex- 
ponent yi is 0.38±0.03. Our value of the critical conductivity 
is in good agreement with the two-terminal conductance, but 
is significantly larger than predicted previously 0]. 

The movement of non-interacting charge carriers in a dis- 
ordered two-dimensional system in the presence of a perpen- 
dicular magnetic field B = (t>h/{ea?) is described by a tight- 
binding model on a square lattice with lattice constant a, and 
0/(27r) is the number of flux quanta per plaquette, 

xy 



xy 



The distribution of spatially Gaussian correlated on-site disor- 
der potentials w„ . with zero mean has been generated as de- 
scribed in Ref . |20| from uncorrected random numbers evenly 
distributed between [—W/V, W/V]. Here, W is the disorder 



strength in units of the transfer term V and 



exp(— [(a: — .t')^ + (2/ — j/')^]/(2k^)) defines the potential's 
correlation length Ip = \/2k which was varied in the range 
< K < /s, where Is — a/^/cfis the magnetic length. 

In contrast to the one-band model llOt 11311 or the network 
model ITsiflTll , the energetical position of the transition points 
are usually not known a priori in the TB model. They depend 
on disorder strength, correlation length, and magnetic field. 
Therefore, careful extensive calculations had to be carried out 
to precisely detect the critical points. A slight imprecision 
influences considerably the results for the conductance and 
the conductivity. We used various methods like calculations 
of the energy and size dependence of level statistics (20], lo- 
calization length iEHi . and Hall conductivity to find the exact 
critical energy. The first quantity was obtained by direct diag- 
onalization whereas the latter two were calculated by means of 
a recursive Green function technique ll22il . The two-terminal 
conductance of a L x L sample with two semi-infinite ideal 
leads attached to opposite sides and periodic boundary con- 
ditions applied in the transverse direction was calculated via 
the Landauer formula g ~ TrtH where t is the transmission 
matrix. For a given disorder realization the conductance was 
numerically evaluated using a well known algorithm I23I1 . En- 
sembles of at least N = 10^ realizations were used {N > lO"' 
in case of L/a = 512) to determine the conductance distribu- 
tion and the corresponding mean critical conductance gc{L). 
The critical conductivity was obtained within linear response 
theory using a recursive Green function method 12^1 . 
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FTG. 1: (color online). The relative deviation of the mean criti- 
cal conductance gc{L) from (7c(oo) vs. reseated inverse system size 
L for a disorder potential with correlation n/a = 1.0, disorder 
strengths W/V = 0.125 (A), 0.25 (O), 0.5 (A), 0.7 (O), 6.0 (V), 
andS = l/(8a^) • h/e. The best fit gives 3c(oo) = 0.6 ± 0.02 and 
a slope y = 0.4. The inset shows the unsealed data for L/a > 40. 



The mean critical conductance gc depends on the system's 
size L, the disorder strength W, and the correlation length 
Ip = \/2k, of the spatially correlated disorder potentials. In 
order to detect the genuine size dependence and to be able 
to correlate our results with those obtained in the presence 
of only one Landau band 

E III E El, 

one has to con- 
fine the parameter range to the situation where the width of 
the disorder broadened Landau bands T{W, k) is small com- 
pared to their energy separation hujc- Stronger disorder en- 
hances those matrix elements that couple neighboring Lan- 
dau bands giving rise to an increasing gc due to extra scatter- 
ing events into the next band, and eventually, shifts the crit- 
ical states to higher energy until they disappear This levita- 
tion of the current carrying states was shown to depend es- 
sentially on the correlation length I21I1 . The size dependence 
of the mean conductance was calculated for various disorder 
strengths W and correlation parameters n/a = 0.3, 1.0, and 
1.5. For disorder potentials with k/u ~ 1.0, which corre- 
sponds to a correlation length of the order of the magnetic 
length {Ip/ls = 1.2533), the rescaled mean conductances are 
shown in Fig.^for system sizes 24 < L/a < 512 and disor- 
der strength 0.125 < W/V < 0.7 and W/V = 6.0. In ad- 
dition, the average over 1610 realizations for systems of size 
L = 1024 and W = 0.5 is also included. The inset shows 
the original data for L > 40 a. The best fit, leading to a mean 
critical conductance gdoo) = (0.60 ± 0.02) e'^/h and an ir- 
relevant scaling exponent y = 0.4 ± 0.02 is also shown in the 
inset of Fig.^Csolid line). Our mean critical conductance is 
in accordance with the results reported for the network model 

(ilEl. 

The rescaling of the data according to gdoo) — gdL) cx 
{C,{W)/L)y is in principle possible only for small disorder 
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FIG. 2: (color online). The critical longitudinal conductivity vs. 
inverse system width M for a correlated disorder potential with 
W = 0.5, K/a = 0.3, and magnetic flux density B — l/(8a^) -h/e. 
The fit to Eq. gives crS^(oo) = (0.58 ± 0.02) e^/h and yi = 
0.38 ±0.03. 



strength. For stronger disorder (1.0 < W/V < 5.0), the true 
scaHng regime cannot be reached for system sizes L/a < 192 
because Landau band mixing, which is obvious from the 
strong overlap seen in the density of states g, causes scattering 
into the critical states in higher Landau bands. This disorder 
induced Landau band coupling may be the reason for the dif- 
fering resuh .g^ = (0.506 ± 0.01)eV/iof Ref.O Increasing 
the disorder strength further, one reaches the situation where 
critical states in higher Landau bands get annihilated by their 
corresponding anti-Chern states m. ForW/V = 6.0, when 
only the critical state of the lowest Landau band survives, scat- 
tering into higher conducting states is not possible so that the 
conductance data again show scaling and match up with those 
obtained for small disorder. A similar behavior is observed 
for weak correlated disorder potentials with n/a = 0.3, where 
Ip = 0.424 a is less than half the magnetic length. In this case, 
we found gdoo) = 0.59 ± 0.03 and y = 0.42 ± 0.03 from the 
rescaled data (not shown). 

Having established the critical conductance in the lowest 
Landau band, we turn now to our results obtained from the in- 
vestigation of the critical Kubo conductivity for disorder po- 
tentials with k/u = 0.3. In Fig.|2lthe dependence of a^^ on 
the width of the system is shown in the range 40 < M/a < 
512. We find in the large size Umit cr^^, = (0.58 ± 0.03) e'^/h 
which is significantly larger than proposed previously. 

The reason for a critical conductance larger than 0.5 e^/h 
and, therefore, for the discrepancy with the suggested be- 
havior involving quantum percolation at the critical point in 
the integer QHE 0, 0, E3l may be associated with the mul- 
ti fractal! ty of the critical quantum Hall eigenstates 
l27l 1281 l29r . a property that was not fully appreciated at that 
time. The diffusion coefficient D, determining the conductiv- 
ity (Jxx = e^gD, is essentially influenced by the multifractal 
correlations of the wave function amplitudes 
Fal'ko and Efetov have shown 1321 13311 that for length scales 




FIG. 3: (color online). The critical conductance distribution func- 
tions P(ln(gf)) from the lowest (O, L/a = 192, W/V = 
6.0, fc/a = 1.0, B = l/(8a^) • h/e) and the second lowest (o, 
L/a = 512, W/V = 0.5, n/a = 1.5, B = l/{32a^) ■ h/e) Landau 
band (LB). The inset shows the relative deviation of the mean critical 
conductance gc{L) for the 2nd LB vs. rescaled inverse system size 
L for a disorder potential with correlation parameter n/a — 1.5 and 
disorder strengths W/V = 0.125 (+), 0.18 (O), 0.25 (V), and 0.5 
(□). The best fit gives (?c(oo) — (0.61 ± 0.03) e'^ /h and a slope 
y = 0.56 ±0.05 



shorter that the localization length, these correlations are also 
present away from the Anderson transition. Using a special 
version of the supersymmetric a model, they provided a for- 
mula relating the diffusion coefficient to the generalized mul- 
tifractal dimensions d{q), 



d{q) 



q 



in^gDh 



(4) 



This relation for unitary symmetry was obtained in the leading 
order in l/{2TigDK) and should hold for q < 2TTgDh. We 
assume Eq. (|4} to be exact near the quantum Hall critical point 
for electronic states with localization length larger than the 
system size ||l33j]. Taking g = 1, we find very good agreement 
with our numerical results presented above. 
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e2 e2 
27r(2 - d(l)) ~ ■ X' 



(5) 



using the information dimension d{l) ~ a{l) ~ 1.739 ± 
0.002 derived H from the precise parabolic f{a{q)) singu- 
larity strength distribution. Over the years, convincing nu- 
merical evidence for a universal f{a{q)) = 2 — {a{q) — 
Q!o)^/[4(ao — 2)] has been accumulated for different quan- 
tum Hall models 111 IH H 13 ES- Within error bars, 
the whole distribution appears to be independent of disorder 
strength, potential correlation length, and magnetic field. It is 
determined only by a sing le value aa = 2.262 ± 0.003 jH. 
Therefore, also d{l) must be universal and so should the crit- 
ical conductivity, if Eq. (|5j was exact. 

We also found the multifractal properties of the wavefunc- 
tions to be independent of the Landau band index for cor- 
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related disorder in agreement with previous investigations 
jyk . Therefore, to check the universaUty of Eq. (|5}, we 
calculated the critical conductance in the 2nd Landau band 
within the same model, but for weaker magnetic field B = 
l/(32a^) • h/e and stronger disorder correlations k/u = 1.5 
(Ip/lB = 0.94). We found ,9c (oo) = (0.61 ± 0.03) eV/i 
(see inset of Fig.|3} from the scaling of the mean critical two- 
terminal conductance, which is in good agreement with our 
results for the lowest Landau band. The critical conductance 
distribution functions P(gf (L)) for the 1st and 2nd Landau 
band of our TB model (shown in Fig.O are almost identical 
and look the same as the one reported for the CC model . 
Conductance values larger that e^/h reflect the fact that more 
than one transmission channels are effective, due to the two- 
dimensional contacts. This is contrary to investigations of the 
critical point-contact conductance f39, '4^ which show a flat 
symmetrical distribution between zero and e'^ /h and, there- 
fore, a mean conductance gc = 0.5 e'^/h. 

A universal critical conductivity a^^. ~ 0.6e^/h should 
be easily observable in experiments on Corbino disks in the 
low temperature limit. Here, the disagreement with the value 
1/2 e^/h should be clearly discernible and, for the first time, 
reveal the multifractality of the critical eigenstates in a trans- 
port measurement. In case of the usual Hall bars, using our 
new value for the critical conductivity, the critical resistance 
Pxx '^xx/ii'^xx)^ + i'^xyf] turns out to be only a Uttle 
smaller, 0.98h/e^ instead of p"^^ ~ 1/i/e^, for the tran- 
sition point in the lowest Landau band, but with 0.23ft,/e^ 
somewhat larger than p'^^ ~ 0.2 h/e^ at the second transi- 
tion. As a further test of the suggested general influence of 
the eigenfunction correlations on the critical conductance, we 
investigated also a two-dimensional disordered system with 
spin-orbit interactions. Using the Ando model, we obtained 
from transfer matrix calculations a critical conductance = 
(1.43 ± 0.05) e^/h and, via direct diagonalization and a mul- 
tifractal analysis of the critical eigenfunctions, an information 
dimension ^"(1) = 1.89. We find again that both quantities 
are perfectly related by a symplectic version of Eq. @ | 32il- 

In conclusion, we have shown that in the critical quantum 
Hall regime both the averaged conductance of square samples 
and the conductivity of very long strips converge to the same 
universal value ~ 0.6e^/h in the limit of large system size. 
This value does neither depend on the considered disorder re- 
alizations nor on the Landau level index, and it is consistent 
with theoretical predictions that take in to account multifractal 
eigenstate correlations at quantum critical points. 
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